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Formation of stable domains filled with strongly correlated coherent quark matter is discussed in general 
terms and is exemplified further in the framework of the Generalised Nambu-Jona-Lasinio model. It is argued 
that such domains, if exist in the Universe, appear dark to an external observer. 
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1. INTRODUCTION 
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^p H . The QCD vacuum is known to be a highly nontrivial 
Oh medium consisting of fluctuations of gluonic fields, vari- 
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>ous condensates, and so on. One of the best studied and 
microscopically understood phenomena related to these 
nontrivial properties of the vacuum is the spontaneous 
breaking of chiral symmetry — a global symmetry in- 
herent to the QCD Hamiltonian, which, however, is not 
respected by its lowest eigenstate, that is by the vac- 
uum. As a result, the entire tower of hadronic states 
built on top of the vacuum lacks this kind of symmetry 
as well. While the effect of spontaneous chiral symmetry 
breaking is expected to play no role for highly excited 
states in the spectrum of hadrons — the so-called ef- 
fective chiral restoration for highly excited hadrons (see 
review |TJ and references therein) — the spectrum of 
low-lying states is strongly affected by chiral symmetry 
breaking. The most prominent example of the manifes- 
tation of chiral symmetry breaking in the spectrum of 
hadrons is provided by the pion, which plays the role 
of the pseudo-Goldstone boson associated with chiral 
symmetry breaking. Another manifestation of chiral 
symmetry breaking is given by the constituent mass, 
of the order of a few hundred MeV, which a light quark 
acquires in the course of its propagation through the 
vacuum. 

The phenomenon of chiral symmetry breaking is a 
highly nonlinear effect, which is described in terms of 
nonlinear gap equations, often referred to as mass-gap 
equations. Such equations are known to possess at least 
two solutions — the trivial solution for the chirally sym- 



metric vacuum |0)o and a nontrivial solution describ- 
ing the physical chirally nonsymmetric vacuum |0). In 
this Letter, we discuss in general terms a physical pic- 
ture which can arise from a possible existence of further 
nontrivial solutions to the mass-gap equation — the so- 
called replicas. The existence of such solutions was em- 
phasised many times in the framework of the Gener- 
alised Nambu-Jona-Lasinio model (GNJL). In particu- 
lar, this was done in Ref. [2] for the quadratic confining 
potential, in Ref. [3J 2] for the linear potential, and in 
Ref. [5] for an arbitrary power-like confining potential. 
In Ref. [3] , similar solutions were considered in different 
approaches, while general arguments in favour of the ex- 
istence of replicas were given in a recent paper [7]. For 
the sake of simplicity, we will assume throughout the 
present paper that only one such excited solution to the 
mass-gap equation exists, which will be denoted as \R). 
Once this is the case, the corresponding state possesses 
the energy density higher than that of the vacuum |0). 
However, since chiral symmetry is broken in the state 
\R), the latter is still energetically more preferable than 
the trivial, chirally-symmetric, vacuum |0)o- Below, it 
is assumed that such a replica state fills a finite-size do- 
main in the Universe, and the problem of stability of 
such a domain is addressed. We start with a qualita- 
tive model-independent consideration, exemplifying it 
further by explicit calculations in the GNJL framework. 

A replica-filled domain at issue is described by just 
two parameters — its radius R (the volume being de- 
noted as V) and the momentum scale A x determining 
the characteristic distance R x at which chiral symmetry 
breaking occurs: R x ~ 1/A X ~ 1 1.5 fm. Indeed, it is 
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legitimate to define A x in this way because the impact of 
chiral symmetry breaking is significantly weaker inside 
the domain than in the outer space, which is filled with 
the unexcited vacuum. For this reason, we can identify 
A x with a typical momentum scale of chiral symmetry 
breaking in the true QCD vacuum. As such, its value 
can be estimated as 250 4- 300 MeV. Accordingly, the 
replica energy density relative to that of the unexcited 
vacuum appears as e = (A x ) 4 > 0. 

Below we impose constraints on the values of _R, 
which ensure stability of the replica-filled domain. 



2. A MICROSCOPIC CONTENT OF THE 
EXCITED VACUUM STATE 

In the true QCD vacuum, spontaneous breaking of 
chiral symmetry leads to the formation of the chiral con- 
densate — a coherent-like state of strongly correlated 
3 Po quark-antiquark pairs. Being yet another genuine 
solution to the mass-gap equation, the excited state \R) 
admits a similar microscopic interpretation. However, 
the way quark-antiquark pairs there get correlated to 
form the condensate is different — while the quantum 
numbers 3 Pq are obviously the same, in order to ensure 
the state \R) to be a scalar, the radial wave functions 
of such pairs appear to be excited. Clearly, while in 
the infinite- volume limit, the |0) and the \R) states are 
strictly orthogonal, for finite volumes V their overlap is 
nonvanishing and reads: 

(0\R) =e- A i v . (1) 

Below, by means of the explicit expressions found in the 
GNJL model, we represent Eq. (TTJ) in terms of the chiral 
angle characterising a given replica \R). 

Thus, we conclude that, from the viewpoint of an ex- 
ternal observer, the replica-filled domain is nothing but 
a localised in space coherent cloud of strongly correlated 
quark-antiquark pairs. 



3. DECAYS OF THE REPLICA STATE AND 
THE MINIMAL SIZE OF THE DOMAIN 

We are now in a position to estimate the minimal 
possible size of a stable replica-filled domain. Indeed, 
Eq. ([1} suggests that the overlap between the replica 
and the vacuum increases as the volume of the domain 
decreases. Accordingly, the lifetime of the domain de- 
creases fast with the decrease of its volume. To quantify 
this effect, one can use the Schrodingcr equation for the 
domain of volume V, 

iSL\R)=U\R), (2) 



where T~L is the QCD Hamiltonian. Namely, we notice 
that: (i) if the excited state \R) decays into the unex- 
cited one, |0), during the time r, the derivative on the 
left-hand side of Eq. can be approximately substi- 
tuted by a finite difference, d\R) jdt -> (|0) - \R))/t; (ii) 
for a sufficiently large volume V , \R) is an approximate 
eigenstate of the Hamiltonian with the eigenenergy Ve. 
Temperature-dependent hadronic contributions to the 
energy density [8, can be shown to be negligible as com- 
pared to e, up to temperatures ~ 25 MeV [7]. 

Finally, projecting Eq. ([5} onto the vacuum state |0), 
using Eq. (TTJ) and the fact that this state is normalised 
as (0 1 0) = 1, one arrives at the following estimate: 
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Thus, we find that the lifetime of the replica-filled do- 
main can be estimated as 
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The condition of stability is obviously formulated as the 
requirement that r is larger than the age of the Universe, 
r > Tu, that is 



e A * y 



(5) 



which yields an estimate for the volume of the domain 
in terms of the dimensionless parameter A x Tjj. 

For A x Tu ^> 1, relation ([5]) can be approximately 
resolved with respect to the volume as 



K>A- 3 ln[A x T i7 ln(A x r [/ )], 



(G) 



which gives for the minimal radius of the domain the 
value 

R min = A" 1 (In [A X T V \n(A x Tu)}) 1/3 (7) 
or, numerically, 



Rn 



5 fm. 



(8) 



Notice that, because of the weak logarithmic, further 
weakened by the cubic root, dependence of the result ([7]) 
on A x Tu, the estimate (JSJ) is very stable against vari- 
ations of the values of both A x and T\j. The obtained 
estimate for i? m i n is of the order of a typical nucleus 
size. 



4. THE MAXIMUM SIZE OF THE DOMAIN 
COMPATIBLE WITH ITS GRAVITATIONAL 
STABILITY 

In this section, we estimate the maximal possible 
size of the replica-filled domain. The existence of such 
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an upper bound for the domain size is a necessary condi- 
tion for the stability of the domain against the gravita- 
tional collapse [7;. Indeed, for a constant energy density 
e at issue, the mass of a replica-filled domain increases 
linearly with its volume. Therefore, the Schwarzschild 
radius r g — 2G ■ (Ve), where G is the gravitational con- 
stant, increases with the size of the domain as 

r g (xR 3 . (9) 

The gravitational stability condition requires that r g < 
R, which leads to an estimate 

i?max ~ (10) 

For a more rigorous analysis, one can apply the Tolman- 
Oppcnhcimcr-Volkoff (TOV) equation for the pressure 
inside a spherically-symmetric domain [7J, and use its 
well-known solution for the case of a constant energy 
density [9 . It yields a more restrictive upper bound 

R max = l/^GA^ ~ 10 km. (11) 

Notice that the obtained estimate for the maximum ra- 
dius is of the order of a typical radius of a neutron star, 
which is also about 10 km. 

A comment is in order here on a possible back reac- 
tion of the gravitational field of the domain on the equa- 
tion of state (EoS) e = (A x ) 4 . Indeed, solving the TOV 
equation together with the EoS, one finds a solution 
which defines unambiguously the gravitational metric 
inside the domain (cf. Ref. [7]). The chiral condensate 
in such a curved space differs, generally speaking, from 
its flat-space counterpart, whose value was initially used 
in the EoS. Correcting the EoS accordingly and plugging 
it again into the TOV equation, one arrives at an iter- 
ative procedure. Fortunately, it turns out that, owing 
to the numerical smallness of the corresponding gravita- 
tional correction to the chiral condensate, the iterations 
terminate already at the zeroth one, that is the EoS 
and the mass-gap equation decouple from each other. 
Namely, the corrected chiral condensate has the form 

(W) = {iji,) [l + O((R x /R) 2 )]. (12) 

This result shows that, since chiral symmetry breaking 
and gravity are effective at different length scales, the 
equations describing them can be clearly separated. In 
particular, one can see that, even for a domain of the 
smallest possible size ([7|). the impact of its gravitational 
field on the chiral condensate is negligible. 



5. REPLICAS WITHIN THE GNJL 

One of the most successful models in studies of var- 
ious phenomena related to chiral symmetry breaking 
is the famous Nambu-Jona-Lasinio (NJL) model |10) . 
This model still possesses a few shortcomings, among 
which the most essential ones are the absence of an in- 
trinsic scale, that shows up in the ultraviolet regulari- 
sation of the loop integrals, and the absence of confine- 
ment. A generalisation of the NJL model was suggested 
in Refs. [5J [TTJ [T^] with the purpose to include con- 
finement, that resulted in the proposition of the GNJL. 
Finally, the direct connection between the 3 Pq quark 
content of the vacuum — in terms of the chiral angle 
tp p (see below) — and the mechanism for the sponta- 
neous breaking of chiral symmetry was shown explicitly 
in Ref. H2]. 

The Hamiltonian of the GNJL reads: 

E = /^,t) Hr v + m )^t) 

+ 1 J d 3 xd 3 y J«(x,t)K^(x - y)JS(y,t), (13) 

where J£(x,t) = , ip(x,t)j l j,^-i^(x,t), and ij){x,t) is 
the quark field. Thus, the interquark interaction is 
parametrised by the quark kernel which contains 
confinement explicitly, in the form of a rising potential. 
The latter can be supplied with various extra terms, for 
instance the perturbative Coulombic one. The confin- 
ing interaction between quarks brings a nonperturbativc 
scale into the problem and models radial excitations of 
3 Po pairs mentioned above, in Sect. 2. Chiral symme- 
try breaking can then be described by the summation 
of loop diagrams for valence quarks (which leads to the 
mass-gap equation), while mesons are obtained from 
the Bethe-Salpeter equation for the quark-antiquark 
bound states. Besides that, GNJL is known to ful- 
fil the low-energy theorems of Gell-Mann, Oakes, and 
Renner (T3], Goldberger and Treiman [13], Adler self- 
consistency zero [15] . the Weinberg theorem [16], and 
several others. We would like to emphasise a universal 
nature of the above low-energy theorems, regardless of a 
particular form of gluonic interactions which result in a 
chiral-symmetry-breaking quark kernel. Consequently, 
irrespectively of a particular form of the confining ker- 
nel K^, GNJL provides a reliable phenomenological 
approach to chiral symmetry breaking, which incorpo- 
rates effects of confinement. Thus, owing to the fact 
that neither qualitative nor quantitative predictions of 
the GNJL depend on a particular form of K ^ , it suffices 
to impose this kernel to be confining and to introduce 
the scale of chiral symmetry breaking, A x . 
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Furthermore, all possible quark-quark interactions 
described by the Hamiltonian (fl3l) include quark self- 
interactions. These self-interactions can be eliminated 
by the use of an appropriate Bogoliubov-Valatin trans- 
formation from bare quarks to the dressed ones. Such 
a transformation can be conveniently parametrised by 
means of the so-called chiral angle <p p (p being the 
relative momentum of a dressed quark-antiquark pair) 

EHH2]: 



E 

p,s 



u(p) = |y i + sin fp 



{otp)\jT 



«o(p), 



v(-p) = -j= + {ap)y/l - sin<^ p J v (-p), 

where a is the colour index of Nc colours. It is con- 
venient to define the chiral angle varying in the range 
—tt/2 < tpp ^ 7r/2 and respecting the boundary condi- 
tions: 



(p p {p = 0) = tt/2, <p p (\p\ ->oo)^Q. 



(14) 



The normal-ordered Hamiltonian ([15} takes the 
form: 

H = E mc + : II, : •://; :. (15) 

and the usual procedure to minimise the vacuum energy 
E v&c is to set the quadratic part : H% '■ to be diagonal. 
Then the corresponding mass-gap equation, 

SE mB [ip]/S<p p = 0, E V3JC [<p] = (0[<p}\H\Q[<p}), (16) 

ensures the anomalous Bogoliubov terms b'd} and db 
to be absent in : H2 :. The mass-gap equation (fTB]) is 
known to possess one trivial solution <p p = and mul- 
tiple nontrivial solutions [2J UJ H3 IH1 H] ■ For the sake of 
simplicity, let us assume the existence of just two non- 
trivial solutions — one describing the vacuum and the 
other defining one replica (more than one replica will 
not change the argument). 

As soon as the mass-gap equation is solved and a 
nontrivial chiral angle is found, the Hamiltonian (fT"5)) 
takes a diagonal form, 



H = E V 



E 



ps ^psi 



(17) 



where E p is the dressed-quark dispersive law, and the 
ellipsis stands for the omitted terms responsible for the 
formation of bound states of quarks. A Fock space can 
be built on top of the nontrivial vacuum |0) upon the 
action of the quark creation operators. 



As it always happens after a Bogoliubov-Valatin 
transformation, the new vacuum contains an infinite set 
of strongly correlated 3 Pq quark-antiquark pairs |12j . 
The same holds true for the replica \R), so it is not 
hard to find explicitly the operator which defines a pseu- 
dounitary transformation from |0) to \R): 



\R) 



|o>, 



(18) 



where C p = b p l[(a-p)ia2] S s'd p l, , with cr's being the 2x2 
Pauli matrices, and Aip p being the difference between 
the vacuum and the replica chiral angles. The operator 
C p creates a 3 Pq quark-antiquark pair with zero total 
momentum and the relative three-momentum 2p, while 
the operator exp[Q — Q*] creates a strongly correlated 
cloud of such pairs. 

Using the commutation relations for the quark op- 
erators b and d 7 one can readily arrive at Eq. ([1} in the 
form: 



(0\R) = exp 



V 



d 3 P 



In f cos 2 



(19) 



The function Aip p is predominantly nonvanishing at 
\p\ < A x , while decreasing fast at larger relative mo- 
menta. For this reason, the integral in the exponent of 
Eq. ((T9)) is C(A X ), and only an overall numerical coeffi- 
cient of the order of unity depends on the details of the 
interquark interaction. 

This way, the scale of chiral symmetry breaking 
A x , initially introduced through the quark kernel K^ v , 
appears naturally in the second-quantised formalism. 
Thus we have obtained a particular realisation of the 
physical picture which was described in the previous 
sections in general terms. 



6. DISCUSSION 

In this Letter, we have discussed some physical con- 
sequences of having more than one realisation of chi- 
ral symmetry breaking in QCD. Indeed, given a highly 
nonlinear nature of the phenomenon of chiral symmetry 
breaking, the possibility for the mass-gap equation to 
possess more than one nontrivial solution does not look 
unnatural — on the contrary, one would need to spe- 
cially arrange for a nonlinear equation to possess just 
one solution. In the literature, the presence of mul- 
tiple solutions to the mass-gap equation was verified 
in a number of NJL-type models, and the correspond- 
ing excited vacuum states were constructed explicitly. 
Among such states, only the lowest one can be associ- 
ated with the genuine vacuum of the theory, while the 



Dark quark domains 



5 



others should be interpreted as some scalar excitations 
on top of it — we use the term replicas for them. In par- 
ticular, by putting forward the conjecture that replicas 
may fill domains in the Universe, we have further ex- 
plored stability conditions for the size of such domains. 
It turns out that the allowed domain radius can vary 
from just a few fermi up to a few kilometres. 

The only way for the domain to be detected by 
an external observer is through its decay to the QCD 
vacuum, with the release of the vacuum energy in the 
form of a cloud of hadrons (predominantly pions), which 
would further annihilate, producing light. Notice how- 
ever that, due to the very strong correlations between 
quarks inside the domain, which result in the exponen- 
tially small overlap ((TJ of the replica and the vacuum 
states, such decays occur quite seldom. And indeed, ac- 
cording to Eq. ([3]), the probability for a macroscopically 
large domain to decay per unit time is exponentially 
suppressed. One is therefore led to conclude that, if 
such encapsulated domains appeared in the Universe at 
its early stages, they would have a chance to survive till 
the present time, remaining however dark to external 
observers. 
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